Scaling relation for the superfluid density in cuprates - origins and limits 
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A universal scaling relation, ps oc o-(Tc) x Tc has been reported by Homes et al. (Nature (London) 
430, 539 (2004)) where ps is the superfluid density and cr(r) is the DC conductivity. The relation 
was shown to apply to both c-axis and in-plane dynamics for high- Tc superconductors as well as to 
the more conventional superconductors Nb and Pb, suggesting common physics in these systems. We 
show quantitatively that the scaling behavior has several possible origins including, marginal Fermi- 
liquid behavior, Josephson coupling, dirty-limit superconductivity and unitary impurity scattering 
for a d-wave order parameter. However, the relation breaks down seriously in overdoped cuprates, 
and possibly even at lower doping. 
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The absence of a recognized theory of high-Tc su- 
perconductors (HTS) has led to a search for universal 
relationships that might provide a guide to the essen- 
tial physics for HTS. The Uemura relationP, Q is one 
such scaling relation, namely Ps(0) oc Tc, where Tc is 
the superconducting (SC) transition temperature and 
Ps(0) = = ^Qe^ns/m* is the superfluid density. Here 
Ug is the density of SC electrons, m* is their effective mass 
and Xab is the in-plane London penetration depth. The 
Uemura relation has been invoked in support of Bose- 
Einstcin condensation of real-space pairs and is gen- 
erally discussed as a test of theoretical models 0. How- 
ever, there is increasing recognition that this relation is 
an oversimplification|^Q, and it breaks down on the 
overdoped side of the SC phase curve 0- 

Recently, a new scaling relation, ps oc ct{Tc) x was 
reported by Homes et where a(T) is the DC con- 

ductivity. This was shown to apply over five orders of 
magnitude including both c-axis and in-plane dynam- 
ics for HTS as well as to the conventional superconduc- 
tors Nb and Pb. The authors suggested this relation 
may provide new insights into the origins of SC in HTS 
materials. We examine this relation and show that it 
is a natural consequence of several quite different but 
well-understood mechanisms, including dirty-limit con- 
ventional SC (Pb and Nb), Josephson coupling along the 
c-axis, marginal Fermi-liquid behavior (optimal and un- 
derdoped HTS), in-plane Josephson-coupled granular SC 
(strongly underdoped HTS) and Abrikosov-Gor'kov d- 
wave pair-breaking for non-magnetic scatterers. Since 
we submitted this work Homes et a/.0 have reported 
some overlapping results but here we examine the doping 
and impurity dependence of this relation showing that, 
in overdoped HTS, it breaks down seriously. It is perhaps 
significant that, here, the cuprates progress towards more 
conventional SC and normal state (NS) behavior. 

Because = /ioe^ris/m* and a = ne^r/m* the cor- 



relation is equivalent to a relationship between scattering 
rate 1/r and Tc given by U/t = 2.7±0.5 x ksTc. Here we 
assume that all the spectral weight associated with the 
free carriers condenses into the ^-function at a;=0, i.e. 
Us—n. In fact, for underdoped cuprates some spectral 
weight remains at finite frequency [lOj. However, we as- 
sume that this weight is related to some other excitation, 
different from the free carriers, which does not contribute 
to a{Lo—Q)^ Ps or the low-T specific heat. 

In considering possible origins for this relation between 
Ti/t and ksTc we note that it is a direct prediction from 
Marginal- Fermi- liquid theory, where [lllll^ 

Ti/t -KksT + huj. (1) 

At low frequency, the scaling relation is almost exactly 
recovered provided r is evaluated at Tc, as required. 

That Nb and Pb fit the scaling line is surprising, but 
it is readily shown that for classical SC in the dirty limit 
Ps is, again, proportional to cr(Tc) x Tc. For a mean-free 
path, ^, and a BCS coherence length, — TT'Vf / {ttA.), 
the effective penetration depth is given by 13] 

Ac// = Al(1 +eo/^)'/' ~ \L{Co/i)'^'. (2) 

Taking £ = vpT and A = l.TGfcsTc, we find that ps is a 
factor of two higher than the scaling line {Ti/t = 5.4 x 
ksTc). In the clean limit there will be large deviations 
below the scaling line because ps remains constant while 
cr(Tc) X Tc can be extremely large. Using data for Pb 
alloys and for Nb alloys we find that, for the Pb 
sample on the Homes scaling line, ^o/£ = 1.2 while the 
two Nb samples have £,o/£ = 0.6 and 2.1. From eq. (2) 
it can be seen that in the intermediate situation where 
^o/^ ~ 1 the above factor of two is compensated. The 
apparent scaling here is understandable but fortuitous. 

Homes et al^ consider Josephson coupling along the c- 
axis and show the proportionality ps oc a xTc but do not 
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FIG. 1: Superfluid density, ps versus cr(Tc) x Tc for 
(a) Yi_a:Ca^Ba2Cu307-5, (b) Bi2Sr2CaCu2 0g+5 and (c) 
Tl2Ba2Cu06+« (overdoped only) and La2-a;Sr2,Cu04. Ar- 
rows indicate optimal doping. In (a) the effects of Zn-induced 
pairbrcaking are shown (crosses) along with the calculated A- 
G pairbreaking effects for constant a (dashed curve) and for a 
reduced by impurity scattering (dash-dot curve). The Homes' 
data lie within ^ ±25% of the scaling line. 



evaluate the proportionality constant. The c-axis pene- 
tration depth, Ac is expressed in terms of the Josephson 
current density, Jc as — fi/ {ifiQdeJc) where d is the 
spacing between SC planes and, for low temperatures, 
Jc = 7rA/(2ei?„)[ll|. Rn = d/uc is the NS tunnelling re- 
sistance between layers. For a d-wave gap we may assume 
that the effective gap parameter, Ae// = Ao/a/2 where 
Aq (~ iAksTc) is the maximum d-wave gap near (tt, 0). 
Thus A-2 ^ 6.79(7rV?ic2)fcsTcCrc giving the relation 
Ti/t = 2.23 X ksTc- If instead we consider a simple rectan- 
gular I-V Josephson characteristic, then 1^ = 2Ao/ei?„ 
and the relation becomes Ti/t — 2.83 x ksT^, again in 
excellent agreement with the scaling curve. 

Thus, whether we consider Josephson coupling, dirty- 
limit SC or marginal-Fermi-liquid SC we recover the ob- 
served scaling behavior. Now we turn to the experimental 
data for the doping dependence of the in-plane dynamics. 

In Fig. 1 the filled squares and diamonds 
show values of Ps(0) plotted versus cr(Tc) x Tc for 
four different HTS materials. The samples, tech- 



nique for measuring Ps(0) and sample-type are as 
follows: (a) Yi_a;Ca2;Ba2Cu307_5 (transverse-field 
muon spin relaxation - /iSR, polycrystalline), (b) 
Bi2Sr2CaCu20g+5 (field-dependent specific heat, poly- 
crystalline), and (c) Tl2Ba2Cu06+5 (^SR, polycrys- 
talline) with La2-a;Srj;Cu04 (ac-susceptibility, oriented 
powders). The /9s(0) data for Tl2Ba2Cu06+5 are from 
ref. 13 and the remainder are from references in 
The techniques noted are well established except 
for the specific-heat method. But /Os(0) has also re- 
cently been measured for Bi2Sr2CaCu208+5 using ac- 
susceptibility 0|. Apart from the heavily underdoped 
region, the data are in excellent agreement with those 
used here from specific-heat. 

The in-plane cr(T) data in Fig. 1 are from DC trans- 
port measurements and for these it is essential to use 
single-crystal measurements. These are available for all 
cited samples except Y^-^,Ca^^a,^?.Gn?. 0^-fi for which 
we have used thin-film data from Fisher llSl. supported 
by more recent data from Naqib et aL[19j. Sources 
for the remaining samples are: Bi2Sr2CaCu208-i-5 from 
Watanabe et aL|20j; for Tl2Ba2Cu06+5 from Tyler and 
coworkers ^ and for La2_2;Sr3;Cu04 from Ando and 
coworkers 23- In all cases we evaluated (t{Tc) at the 
reported doping states, fitted to an appropriate func- 
tion of doping and interpolated to determine cr(Tc) at 
the doping state relevant to the available ps data. For 
Bi2Sr2CaCu208+5, the (j{Tc) data had to be extrapo- 
lated to higher doping levels. In Fig. 1 optimal doping 
(maximal Tc) is indicated by the arrows and critical dop- 
ing, where the pseudogap closes, is the point of maximum 
Ps[3- The overdoped region lies further to the right. 

Our main conclusion is that each compound exhibits 
a large departure from the "universal" scaling line in 
the overdoped region. In the optimal and underdoped 
regions, where we can directly compare with Homes 
et al, the situation is less clear. In some cases our 
data differs markedly from theirs which severely over- 
estimates Ps for Bi2Sr2CaCu208 relative to data from 
specific heat, ac susceptibility and fiSK and underesti- 
mates Ps for Tl2Ba2Cu06. For this compound, and for 
La2_2;Sr2;Cu04, their a{Tc) value is about half that ob- 
served from DC transport in single crystals 21). The 
source of these discrepancies may possibly be attributable 
to the extrapolation to a; = and a; = 00 required for 
the Kramers-Kronig transformation of the optical data. 

Bearing in mind these differences our underdoped re- 
sults show a consistent trend. For Bi2Sr2CaCu208, 
Tl2Ba2Cu06 and La2-£i;Sra;Cu04 the scaling line is not 
reached until near, or below, 1/8*'' doping. On the other 
hand Yi_2;Caa;Ba2Cu307, though not hnear, roughly fol- 
lows the trend of the scaling line across the entire under- 
doped region. However, there are two reasons why this 
could misrepresent the situation. Firstly, compared with 
other bilayer cuprates with similar Tc this compound gen- 
erally exhibits a higher superfluid density arising from 
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the contribution of the chain layer. This may be approx- 
imately balanced in the scaling plot by the additional 
conductivity also arising from the chains and so we put 
this effect to one side. But secondly, the conductivity 
here is for thin films and thus is probably diminished rel- 
ative to single-crystal data. All the data points may lie 
further to the right, perhaps by a factor of the order of 
1.5. Then the behavior would be rather similar to that 
shown by Bi2Sr2CaCu208+5 and La2-2:Sr2,Cu04. It may 
be that there is no correlation with the scaling line until 
the lowest doping state around p = 1/8 and lower. 

What can we conclude from this? One implica- 
tion may be that in this heavily underdoped region 
the in-plane physics is rather similar to the out-of- 
plane physics, namely that of Josephson-coupled SC do- 
mains. Here the electronic state is probably granular, 
with SC patches separated by insulating Josephson bar- 
riers. Such a model has been proposed even for opti- 
mally doped cuprates and STM suggests the presence of 
nanoscale inhomoereneitv jj^t .24] . However this has been 
questioned 25], and in particular the electronic state seen 
from the perspective of NMR and specific heat seems 
to be homogeneous for p > 0.12[2^. Below that dop- 
ing state these techniques reveal the presence of normal 
quasiparticlesji^. Here the electronic state definitely ap- 
pears to become spatially heterogeneous. 

We turn now to ask the question as to whether the 
universal scaling relation may also apply in the case of 
impurity scattering. This behavior can be understood in 
terms of generally-accepted theories for the pair-breaking 
effects of Zn impurities, namely the suppression of Tc by 
non-magnetic scattering, which for a d-wave supercon- 
ductor suppresses Tc according to the Abrikosov-Gor'kov 
(A-G) formula and ps according to the work of Maki and 
coworkers 28J. When we add in the effects of scattering 
on a{T) (described by Matthiessen's rule) this roughly 
reproduces the observed scaling behavior for optimally- 
doped samples but, again, we expect large deviations for 
overdoped samples for which conductivity data as a func- 
tion of Zn doping are not yet available. The theory can 
be developed along the following lines. 

For a non-magnetic scatterer d-wave SC is rapidly sup- 
pressed and the effects on Tc, ps and cr as a function of 
impurity concentration, y, are governed by the density 
of states (DOS) at the Fermi level. These effects have 
been shown to be fully consistent with thermodynamic 
dataji^. The reduction in Tc is given by: 



in(rc/Tco) = ^[1/2 + r/rj - v[i/2] 



(3) 



where tpi^] is the digamma function, Tco — Tc{y—0). For 
unitary scattering, F = ni/'Kg{Ep) is the pair-breaking 
scattering rate and Tc is its critical value for which Tc 
is suppressed to zero. Here g{Ep) is the DOS per spin, 
ni{= ay/abc) is the density of impurities, a is the num- 
ber of Cu02 planes per unit cell and a, b and c are 
the lattice parameters. As before^^, we adopt the 



strong-coupling scenario Tc = 1.65kBTco and determine 
g{Ep) from the electronic specific heat coefficient, 7, us- 
ing 7 — {2/3)Tr'^k'^g{Ep). The following linearized form 
of the A-G equation is valid up to about (2/3)rc: 

Tc/Tco = 1 - 0.69F/Fc = 1 - (OMaR/jTco) x y, (4) 

where R — Na^b- In the overdoped region 7 is con- 
stant, independent of temperature and doping, so that 
the slope, dTd dy^ in eqns. (3) or (4) remains constant. 
As a practical measure 29], the value of 7 was taken as 
S/T evaluated at Tc, where S is the electronic entropy. 
This is just the average of 7(T) between T = and Tc- 

The calculation of the reduction in superfluid density 
is more complex and is carried out nume ricallyHI. How- 
ever, we find an excellent approximation to within a few 
percent across the entire pair-breaking range using: 



Ps/PsO 



Tc/T, 



cO 



l + 1.7{l-Tc/Tco)' 



(5) 



This non-linear relation reflects the fact that ps is ini- 
tially suppressed much faster than Tc 5]. Eqns. (3) 
and (5) are universal relations and if impurity scatter- 
ing were to have no effect on a then the effect of Zn 
substitution would be given by the dashed A-G curves 
shown in Fig. 1(a) for the optimal and most overdoped 
Yi_2,Ca2,Ba2Cu307_.5. 

In fact, the resistivity is increased by impurity scatter- 
ing according to a linear Matthiessen law. The resultant 
decrease in (7{T) seems to effectively linearize the simple 
A-G curve in the scaling plot. We illustrate this using the 
experimental data shown by the crosses in Fig. 1. These 
are data for an optimally-doped sample where the reduc- 
tion i n Ps and Tc is from /iSR measurements by Bernhard 
et aL|3l|, and the cr(T) data for Zn substitution at the 
same doping state are from thin-film samples studied by 
Naqib et al 19]. These trend back linearly to the origin 
in reasonably good agreement with the scaling curve. We 
now seek to justify this quantitatively. 

The residual planar resistivity arising from impurity 
scattering is given by|3^ A{h/e'^){ni/n)sin'^So where 60 
is the s-wave scattering phase shift which, in the unitary 
limit, we take to be tt/2. Thus the total resistivity is 
m* / {ne'^r) + 4c{h/ae'^){ni/n) and a becomes 



ne^r/m* 



1 + {ATiT/abm*) x y 



(6) 



where t is evaluated at Tc{y). 

There are two impurity effects to be considered. 
Firstly, there is the Matthiessen term in the denominator 
of eqn.(6) and, secondly, cr(T) must be evaluated at Tc 
which is itself reduced by scattering. This second effect 
is accommodated by the fact that, in eqn. (7) below, 
the only two variables in the coefficients are t and Tc 
and they appear as the product tTc- Thus, provided the 
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material lies on the scaling curve in the absence of impu- 
rities, this product is a constant. Proceeding, the explicit 
variable y in cqn. (6) may be eliminated by substituting 
from eqn. (4) to give: 



a(0)T,o X (T,(y)/T,o) 



(7) 



l + (/3-l)(l-T,(y)/reo)- 

where j3 = Ah'^ToTcQ/O.SGaRabm* and tq = r evaluated 
at Tc{y ~ 0). Taking m* « 2me, as observed for nodal 
quasiparticles^^, we have only to note that the right 
side of eqn. (7) is of the same form as eqn. (5) with the 
coefficient (/3 - 1) = 1.9. Thus 



(^{y)Tc{y)l Ps{y) ~ a{0)Tco/pso, 



(8) 



and the impurity scattering data remain on the scaling 
curve if, in the first instance, it lies on the curve in the 
absence of impurity scattering. This seems to explain 
the crosses shown in Fig. 1. This calculation shows that, 
with underdoping, the curve would steepen due to the 
further opening of the pseudogap (and the associated re- 
duced entropy at Tc) while with overdoping it would flat- 
ten (due to the closing of the pseudogap). It would be 
nice to test this but we do not have combined Ps{y) and 
cr(y) data as a function of Zn concentration for heavily 
overdoped samples. Nevertheless, we have calculated the 
expected effect of Zn substitution using eqn. (7) and this 
is plotted in Fig. 1 by the dotted line. This shows a 
fairly linear behavior but which is displaced far from the 
"universal" scaling curve. 

In conclusion, we have examined the scaling relation, 
Ps (X <j(Tc) X Tc, and find this is equivalent to the relation 
Tl/t = 2.7 ± 0.5 X ksTc which could equally arise from 
(i) conventional dirty-limit SC, (ii) marginal- Fermi- liquid 
behaviour, (iii) Josephson-coupling along the c-axis of 
the Cu02 planes, or (iv) unitary-limit impurity scatter- 
ing. We conclude that the Pb and Nb samples correlate 
because these samples are near the dirty limit. Interest- 
ingly, the Homes scaling line corresponds to the mean 
free path just above Tc being approximately 2 x (us- 
ing the maximum gap « 2.38fcBTc for a weak-coupling 
d-wave SC). Even more intriguingly, in terms of stan- 
dard microscopic BCS theory it corresponds to the 
Gor'kov kernel having a range near Tc equal to the 
mean free path just above Tc. 

When we examine the doping dependence of the in- 
plane Ps and o'(Tc) x Tc we find a total breakdown of 
the scaling relation across the overdoped region. This 
shows that the mean free path at Tc is now much larger 
than both of the above-noted coherence lengths. This 
breakdown may extend into the underdoped region, al- 
lowing for the fact that the only exception is for thin- 
film Yi_a;Caa;Ba2Cu307_5. The recovery of the scal- 
ing behavior at very low doping may arise from simi- 
lar physics to c-axis Josephson coupling, namely, weak- 
linked patches of superconducting domains in the inho- 
mogeneous strongly-underdoped state. 
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